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We derive a functional equation for the mean first-passage time (MFPT) of a generic self-similar 
Markovian continuous process to a target in a one-dimensional domain and obtain its exact solution. 
We show that the obtained expression of the MFPT for continuous processes is actuaUy different 
from the large system size limit of the MFPT for discrete jump processes allowing leapovers. In 
the case considered here, the asymptotic MFPT admits non-vanishing corrections, which we call 
residual MFPT. The case of Levy flights with diverging variance of jump lengths is investigated in 
detail, in particular, with respect to the associated leapover behaviour. We also show numerically 
that our results apply with good accuracy to fractional Brownian motion, despite its non-Markovian 
nature. 

PACS numbers: 05.40.Fb 



The first-passage time (FPT) of a random walker to a target point is an important concept for the modelling and 
interpretation of stochastic processes pQ [5] . This interest in the concept of FPTs is motivated by the crucial role 
played by FPTs in various contexts, including diffusion limited chemical reactions [3HS], the spreading of diseases [5], 
or target search processes [TUlOj. A general formalism for the calculation of the mean first-passage time (MFPT) 
of a scale-invariant random process Vt in a confined domain of volume V has recently been derived in the large V 
limit in fTTJ [T^] , and the full distribution of the FPT has been obtained [13] ■ These results highlight the asymptotic 
dependence of the MFPT (T) and its higher order moments on both the source-to-target distance r and on the 
confinement volume in the case of scale-invariant processes, which can be characterized by the walk dimension dw 
(defined by (r^) cx t^/'^'^) and the fractal dimension d/ of the support of the random process [T4]. In particular, in 
the case of compact (i.e., recurrent) processes, for which > df,\t was shown that 



where only the constant B depends on the process. For a Brownian walker the first passage at a given point x, which 
hereafter denotes the first arrival at exactly x is equivalent to the first crossing of the value x. This is not generally 
true. For instance, in the case of Levy fiights with scale- free distribution of jump lengths [T5til7j . pronounced leapovers 
across the threshold value occur. Thus, the first passage at the point x becomes considerably less likely than the first 
crossing of the value x [TS]. For Markovian Levy flights, the first crossing falls into the class of processes governed by 
the Sparre Andersen universality |19) . 

Discrete time Markovian jump processes are defined by an elementary jump distribution w{v), which, at each 
discrete time step n, renders the probability that the random walker makes a jump of length r, which we here assume 
to be discrete. For a generic distribution w, the searcher is in principle allowed to jump across the target, and perform 
what we define here as a leapover, as opposed to the case of nearest neighbor random walks. In the case of Euclidean 
spaces, symmetric jump processes are subject to the generalized central limit theorem [201 [H], and have a well defined 
continuous limit which is scale-invariant. More explicitly, let -P(r, n) denote the discrete time and space propagator 
of the process in infinite space, starting from r = at time n = 0. Then we know that for the rescaled quantity 
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the convergence to 



limpe(x = re, t 



(3) 



2 



is fulfilled in the continuum limit e — ^ with x and t fixed by virtue of the central limit theorem. Here g is a scaling 
function. If the jump distribution w{t) has a second moment (r^), the limit process is a Brownian motion with = 2 
and g is Gaussian; otherwise the limit process is of Levy stable type with < dyj < 2, and g is a Levy stable law 
of index d^,. In this continuous limit the asymptotic theory applies fllj, for compact processes yielding the exact 
asymptotic form ([T]), where B can be explicitly determined as a function of g. These results show in particular that 
the MFPT in the continuous limit becomes independent of the existence of leapovers. 

In terms of the variables r and n of the discrete jump process, the continuum limit is equivalent to taking both r 
and n large, with r'*"'/n fixed. Hence, for discrete jump processes, Eq. ([ij produces actually only the leading term of 
a large r expansion, and one should write more precisely 

lim ^ = Br'^^-'^f + o{r'^"'-'^f). (4) 

We show in this paper that the subleading term of this expansion actually strongly depends on the existence of 
leapovers, and may yield non-vanishing corrections to the rescaled MFPT. This residual MFPT can have important 
consequences in the context of both search processes and numerical simulations of first-passage times of random walks. 

As a specific important example we discuss extensively the effect of leapovers on the MFPT in the context of 
one-dimensional (II?) Markov jump processes. By means of a functional equation, we first derive an exact expression 
of the MFPT for continuous time and space scale-invariant Markov processes valid for any starting position of the 
walker and confining domain size. Next, using the method of pseudo Green functions, we present an exact expression 
of the MFPT for arbitrary discrete time ID Markov jump processes, which was previously derived in Ref. [1, 22 . We 
show that in the large volume limit, the leading term of the MFPT is the expected continuous limit of Eq. ([T]), and 
we calculate exactly the subleading term of this expansion. We show that this correction is indeed crucial in the case 
of random walks with leapovers, since it does not vanish in the limit L — > cxi and r — > 0. 

We consider a walker performing random independent jumps in a confined ID system of size L (with periodic 
boundary conditions), and address the question of determining its mean first-passage time (Tts) = (T(r, L)) to a 
given target site T as a function of its starting site 5, where r stands for the source-to-target distance \ST\. To 
proceed, we denote by Wji the jump probability from discrete site i to j, where Wji is assumed to be symmetric 
(wji = Wij) and translation- invariant (wji is a function |i — j| only). It is easily seen that the confined problem with 
periodic boundary conditions is equivalent to an infinite line with regularly spaced targets at positions kL with G Z, 
a property which will be used below. 



II. CONTINUOUS TIME AND SPACE SCALE-INVARIANT JUMP PROCESSES 



In this section we derive an exact expression for the MFPT of a ID scale-invariant Markov jump processe in the 
limit of continuous space and time. We start from the definition in discrete time and space given above, and pass to 
the following classical backward equation for the MFPT [5]: 

CO 

J2 (Tn)(«izs - S,s) = As{Tts) = -1, (5) 

i— — oo 

completed by the condition (T(r = kL, L)) — for all k G Z. It is easily seen that the following form of the MFPT 

^^(^^ ^ (T(r, 2L)) + (T(r + L, 2L)) - {T{L, 2L)) 

satisfies both Eq. ([s]) and the boundary conditions above. Note that by definition, {T{r,L)) is L-periodic. This 
provides a functional equation, which can be solved as follows. We assume that the jump process has a well defined, 
self-similar continuous limit, characterized by a walk dimension d^,. We then define the continuous limit of the MFPT 
as 

eiy, I) = hm e"- (T(r = y/e, L = l/e)), (7) 

e->0 

where the starting position y and system size I are fixed. In particular, we assume that the process is compact 
{diu > df ~ 1), such that in the continuous limit the MFPT to a point-like target is finite. In this limit self-similarity 
implies 

e{y,l) = l''-9{x^y/l,l), (8) 
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and therefore 



1) = 2'--' 1) + ei^, 1) - ei^^ 1) 

so that we now have to solve the foUowing functional equation 

fix) - 2-*™-! 



x + l 
2 



We know that /(x) = /(I — a;) by symmetry, and that /(O) = /(I) = 0. Differentiation of this equation yields 



fix) ^2^--^ (^/'(|)+/'(^: 



(9) 



(10) 



(11) 



together with the symmetry condition J'{x) — —/'(I — x). It can be checked directly that equation (11) is solved by 
the Hurwitz zeta function Q{2 — diu,x), the latter being defined by the series 



oo ^ 

Cis,q) = Y] 1—, — 



(12) 



ri=0 



Note that C(2 — d^,, 1 — .t) is also a solution. In order to satisfy the symmetry condition, a solution of Eq. ( 11 ) satisfying 
f'{x) — —/'(I ~ a;) is C,{2 — d^,x) ~ ({2 — dw, 1 — x). A solution of the original equation (10) satisfying /(O) = is 
then 



f(x)= / {C(2-d^,y)-C{2^d^,l-y))dy. 



Furthermore, using the fact that the Hurwitz zeta function satisfies the relation 

d({s,x) 



dx 



-sC(s + l,x), 



one can simplify Eq. ( 13 ) and write 

1 



(C(l - d^, x) + C(l - 1 - x) - C(l - 0) - C(l - 1)) . 



dw ~ 1 

The MFPT eventually becomes 

e{y, I) = AL"- (C(l - |) + C(l - 1 - |) - C(l - rf»,0) - C(l - d^, 1)) 



(13) 



(14) 



(15) 



(16) 



where A is a multiplicative constant which remains undetermined at this stage. For d^ — 2 (Brownian case), 
we retrieve the classical expression, using C,{—n,x) = —Bn+i{x)/{n + 1), where i?„ are the Bernoulli polynomials 
(^2(2;) ^x^ ~x + l/&): 

fix) = c(-i,.T) + c(-i,i-a:)-C(-i,o)-C(-i,i) 

= -^iB2{x) + B-2il - x) - B2{0) - 82(1)) 
= x{l — x). 

Comparison of this expression to the well known result 9{y, I) = y{l — y) yields A = 1 in this Brownian case, with 
D = 1/2. 

In the general case d^ ^ 2, A can be calculated by using the exact asymptotic result in the large / limit of Eq. ([T]) 
(see Ref. [TT] for details). For y ^l, the asymptotic behavior of Eq. (16) readily produces 



lim«=V^-. 



Z— f 00 / 

Comparison with the general result of Ref. [TT] then delivers 

du 



A = 







(17) 



(18) 
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where the infinite space propagator of the process p(r, t) satisfies: 

p{r,t)=t-'^''-g{r/f^''-), 



(19) 



and g*{u) — g{0) — g{u) tends to when u — )• 0. 

Because of the generaUzed central limit theorem invoked in the introduction, a generic continuous time self-similar 
Markov process with 7^ 2 is a Levy process, whose infinite space propagator reads 



p(r,t) 

where c is a scaling parameter. We subsequently deduce 

1 



cos(rx)e' 







°° I -( 
— cos(w)e u 
u 



dv. 



and 



9*iu) 



1 

- (1 - cos{v)) e u dv. 
u 



(20) 



(21) 



(22) 



The constant A can then be explicitly calculated and reads: 

du 



A = 



A = 







1 — cos(w) 
tt(cv)'^^ 



^1"- (l-cos(t;))e-M''""di; 



dv 



2c<i»r(l + (i^)cos 



(2 - dw)T: 



(23) 



Equations ( 16 1 and ( 23 1 provide an explicit and exact expression of the MFPT for all r and L, valid for any continuous 
time and space scale-invariant Markov process. 

In order to check equation (16), we simulated Levy fiights with various values of d^,. A Levy flight is here modeled 

;h step-size 

w{x — n) 



as a discrete-time random walk, where each step-size is a random integer x, distributed according to the distribution 

1 



2C(l + d„)|7i|i+'^" 



(24) 



with w{x = 0) = 0. The random walk takes place on a ID lattice of size iV, in discrete time, and with periodic 
boundary conditions. The target is located at the site of origin = N. This processe converges in the continuous 
limit defined in introduction to a Levy process, whose infinite space propagator is given by Eq. (20 1 with 



2d„ 



ndu 



(25) 



r(d„)C(i + d.,) 



We compared the exact expression (16) obtained in the continuous limit with our discrete simulation. Figs. [T] and [2] 
show that as the lattice size N grows, the simulated discrete MFPT approaches the exact continuous limit, for different 
values of d^ (for each dw, the multiplicative constant A is computed using ([23])). Note the very slow convergence to 
the exact solution in the continuous limit. We show below that leapovers lead to non vanishing corrections of the 
MFPT to this continuous limit in the large system size limit, which we define as the residual MFPT. 



III. DISCRETE TIME AND SPACE JUMP PROCESSES 



To investigate the convergence of the discrete process to its continuous limit, we compute analytically the MFPT for 
a generic discrete time and space Markov jump process, regarding jump processes with and without second moments. 
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FIG. 1: Mean First Passage Time for a Levy flight, for several sizes of the ID lattice (200, 400, 800, 1600, 3200 and 6400 sites, 
from black [bottom] to light grey [top] as the size grows) compared to the theoretical expression ( 16 1 (magenta dotted line, the 
A constant being given by Eq. (|23[)), and for = 1.25. 




FIG. 2: Mean First Passage Time for a Levy flight, for several size of the ID lattice (200, 400, 800, 1600, 3200 and 6400 sites, 
from black [top] to light grey [bottom] as the size grows) compared to the theoretical expression of Eq. ( 16 1 (magenta dotted 
line, the A constant being given by Eq. (23l), and for — 1.5. 



We use the formalism developed by Montroll [55] (see also Hughes [T]). First we define the structure function A(fc), 
the Fourier transform of the jump distribution (24): 



X{k) = ^ i(;(a; = 71)6*"^ 



(26) 



Denoting by Pt(fc) the Fourier transform of the propagator at discrete time t in infinite space, we obtain Pt{k) ~ A(fc)*. 
The infinite space Green's function then reads 



00 CO ^71- 

GJ^ = 5]P(j,t|*,0) = ^- / cosik\i-j\)XikYdk 
t=o t=o •^'^ 



1 rcoa{k\i-j\) 



ttJo 1 - A(fc) 



dk. 



(27) 
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Recalling that there exists a set S of targets located at sites nL {n € Z) it is useful to define 



n— — CX3 

1 ^1 cos(^) 



rn — l V 

It can then be checked directly that the following expression 



27rm ^ 
L ) 



L-1 



(29) 



1 — cos 



(T(r, L)) = L{G^^ G^r) = ^ , _ u^A) (30) 

m=l L ) 

represents an exact solution for the MFPT of a ID discrete time and discrete space symmetric random walk, since it 
satisfies both Eq. ([5| and the boundary condition (T(r = nL, L)) = for all n G Z. In what follows we distinguish 
(i) the continuous limit 0{y,l) of this exact expression, defined by Eq.([7]), and (ii) the scaled MFPT in the large L 



limit, which reads according to Eq. (30) 



(T(r,L)) 1 ri-cos(fcr)„ 
lim ^ \' " EE T r = - / ^jYY^dk. 31 

L^oo L TT Jq I - X{k) 



While it is clear that the continuous limit 9{y, I) must yield the exact result of Eq. ( 16 ) derived independently above, 
we will show below that the scaled MFPT r(r) can differ from this continuous limit as soon as the jump distribution 
allows leapovers. The residual MFPT corresponds to the non- vanishing, subleading term of T(r) in the large r regime, 
which can be obtained by considering the fc — > behavior of A(fc). We consider different examples below, corresponding 
to jump distributions with and without second moment. 



A. Residual MFPT for jump distributions with second moment. 

We first consider a random walk whose jump distribution possesses a second moment. As a concrete example we 
employ the exponential law defined by 

P(„) = M_l)e-/3|"l (32) 
for n > 1 and P(0) = 0, where /3 is a scaling parameter. The second moment of this jump distribution becomes 

(n^) = ± n-P{n) = (33) 

n— — oo 

SO that the continuous limit of this jump process, according to Eq. ([3]), is a Brownian motion with diffusion coefficient 
D = (rt^)/2, and therefore dm = 2. In this continuous limit one obtains straightforwardly the following form of the 
MFPT, 

eiyj) = linr e'^-(T(r = y/e,L^ l/e)) = ifj^— ^y(; - y). (34) 
We now apply the formalism of the previous section and first compute the structure function A(fc), obtaining 
A(fc)= y P(a; = n)e"'" = (e'^ - 1) y cos(nfc)e-'^" = (6-3 - 1) cos(fc) - 1 



n— — OO 
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An exact expression of the MFPT is then given by (30), for which we find 

(T(r,i)> 



L-l 



^ l-cos(^) 



27rm ^ 
L ) 



1 1 



f 2'jTmr \ 
'\ L ) 



^ ' m— 1 



COS 



27rm ^ 



X ((e'' - 1)2 + 26^^(1 - cos (^))) 



1) 



6/3(6/3 + 1) 

+ - 



:'3(e/^ + 1) 



r(L — r) 



i. 



(36) 



While this result has the expected continuous limit (34), it reveals a significant difference from the scaled MFPT, 



containing a non-vanishing residual, subleading correction. More precisely we find that 



lim 



(T(r,L)) 



T(r) 



1) 



2e/^ 



e/3(e/3 + l) e^(e^ + l)' 



(37) 



where we define the non-vanishing subleading term 2e^/(e^(e'^ -I- 1)) as the residual MFPT. 

Actually, the residual MFPT for a generic jump distribution with second moment can be fully derived. We write the 



small k expansion of the structure function as A(fc) 
as above. Then one has 



1 — Dfc^, where the diffusion coefficient is given by D = (ri^) /2 



r(r) 



1 

TT 

r 

2D 







1 — cos (fcr) 
Dk^ 

B + oil), 



dk + B + o{\) 



where the residual MFPT is given by 



B 



l-A(fc) 



dk 



(38) 



(39) 



This reveals the crucial difference between the large system size limit on the one hand and the continuous space and 
time limit for MFPTs of jump processes on the other. Our above analysis (see Eq. ([39|) suggests that a residual 
MFPT exists in the large system size limit as soon as leapovers are allowed. In particular the residual MFPT depends 
on the full jump distribution and not only on its second moment. For example, the case of a nearest neighbor random 
walk, for which X{k) — cos{k), yields a vanishing residual MFPT, as can be checked from Eq. (39). 



B. Residual MFPT for Levy flights 



We now turn to jump distributions with infinite second moment. The continuous limit of such jump processes is 
a Levy flight by virtue of the generalized central limit theorem. We start with the example of discrete space Levy 
flights, defined by the distribution (24 1, for which one has 



A(fc) — w{x — n)e' 



1 

CKTI) 



E°° cos(n/c) 



(40) 



We first note that using this expression in the exact result (30) yields in the continuous limit the result (16) and (23), 



as expected. We consider below the large L limit of the MFPT and will therefore make use of the following small k 
expansion [1] 



A(fc) = 1 - aifc''"' + aa/c^ + ©(fc"*). 



where 



"1 



2C(d™ + l)r(d^ + l)sin(^)' 



(41) 
(42) 
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and 



C(2 - d^)r(2 - du,) ( TTd. 



a2 = n A COS — — . (43) 



We then compute the following large r expansion of the scaled MFPT: 

dfc + 0(r^''™~^). (44) 



1 ri-cos(M.. , 



For any value of du, e ]1, 2[, we obtain the following expression of the leading term 



aiTT 



oo 



r(r) ^ ■ / l^^dfc 



(45) 



2air(l + d^) cos 



7r(2 - d„) 



We therefore obtain exactly the result (!]) derived in the continuous limit, with the multiplicative constant given by 
(23). In particular the leading term of the MFPT does not depend on the exact form of the discrete time propagator, 
but only on its continuous limit. 

To assess the subleading term, we have to distinguish between the cases d^j > 3/2 and d^ < 3/2. For the latter 
case, we obtain a correction term to the scaled MFPT behaving as 0(1), 

T(r) = . , ,s r'--' +B + o(l), (46) 

2air(l + d„)cos(^^^=^j 

where the subleading term B, defining the residual MFPT, is a constant, that does not vanish in the limit L — >■ oo, 

1 f 1 1 1 

B=- rvTT TTT-rffc — ^. (47) 

ttJo l-A(fc) aifc^™ ai(d„ - l)7r'^- ^ ' 

The residual MFPT B depends on the full k expansion of A(A:), and thus on the precise form of the propagator. If 
the leading term is the same for all jump processes sharing the same limiting propagator, the residual term depends 
on the shape of the jump distribution, as was found in the case of jump processes with finite second moment. The 
validity of expansion (46 1 is verified numerically in Fig. [3j 

For the case d^ > 3/2 the third term in the fc expansion of A(fc) has to be taken into account. We now need to 
compute a second subleading term in the large r expansion. 



1 Z"' 1 - cos (fcr) „ \ r \ - cos (fcr) „ , , 

T{r) = - I rj-^dk^- — ^-^dk + B + o{l) (48) 

0L2 



a. 



d^.-l , i'^d.yj — 2)q!2 2d„-3 



2ai7rr(l + d^) cos (li?^) 2a2r(2d,, - 1) cos(^d^) 



Importantly, this expansion yields again a residual term B, which does not vanish for small r, and which depends on 
the full jump distribution: 

R = 1 r f ^ - ^— - "2 ^ AU _ 1 «2 , . 

ttJ^ \l-\{k) a^kd^ ajk^d^-y ai(d„ - l)^'^- aU2d.^ - 3)n^'^^-2' ^ ' 

The validity of expansion ( 49 ) is verified numerically in Fig. |4] 

To emphasize how critically the residual MFPT depends on the jump distribution, we consider the example of 
a Levy flight in discrete time, but now in continuous space. At each time step, the random walker makes a step 
distributed according to a symmetric Levy stable distribution of index a — dyj. The structure function X{k) is then 
directly the characteristic function of the Levy distribution, 

\{k) = exp(-|cA:|'^") = 1 - [ckf- - ^ + 0{k''^-). (51) 
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200 



FIG. 3: Mean First Passage Time for a Levy fl ight on a ID lattice of size 200, with = 1^25, in discrete time (black circles) 
compared to the theoretical expression of Eq. ( 30 1 (red line) , to the rough approximation ( 45 1 (green dotted line) and to the 
approximation (46 1 (blue dashed line). The inset shows the difference between the three theoretical expressions (Tth) and the 
simulation (Tsim). 




200 



FIG. 4: Mean First Passage Time for a Levy flight on a network of size 200, for ~ 1.75, in discrete time (black circles) 
compared to the theoretical expression ( |30[ ) (red line), to the rough approximation of ( |45[ ) (green dotted line), and to the 
approximation (49 1 with the first subleading term (blue dashed line) and with two subleading terms (magenta plain line). The 
inset shows the difference between the three theoretical expressions (Tth) and the simulation (Tsim). 



The leading term of the MFPT is then the continuous Umit obtained in Eqs. (171 and (23), with ai = c However, 



the subleading term in this case is always a constant residual MFPT, whose va ue depends on the specific shape of 
the jump distribution. This means that the scaled MFPT admits difi^erent large r expansions for Riemann walks 
(distribution (24)) and continuous space Levy flights, which depend on the full jump distribution and not only its 
asymptotics. Importantly, we found that in both cases there exists a non-vanishing residual correction to the MFPT. 
These results, together with previous section, suggest that a non- vanishing residual MFPT exists as soon as the jump 
distribution allows leapovers. 
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FIG. 6: Splitting probability for Levy flights, for several size of the discrete network (200, 
from black to light grey as the size grows) compared to the theoretical expression of Eq. 
~ 1.5. 



400, 800, 1600, 3200 and 6400 sites, 
( 53 1 (magenta dotted line) , and for 



IV. SPLITTING PROBABILITIES 



The above formalism can be extended to further first-passage observables, including observables involving several 
targets, following the method developed in Ref. [23J. As an illustrative example we consider the case of splitting 
probabilities. More precisely, we assume that the jump process takes place on a ring of length 2L, with a target Ti 
at position r = and a target T2 at position r = L. Note that in the case of processes with leapovers, first-passage 
problems involving N > 2 targets do not reduce to problems with two targets. We denote by Pi the probability that 
the walker hits Ti before ever reaching T2. It is known that Pi can be exactly expressed in terms of MFPTs as follows 

^ ' ' 2(T(L,2L)) ^ ' 

Indeed, it can be easily seen that this form satisfies the boundary conditions Pi{0,L) — 1 and Pi{L,L) = 0, as 



well as the backward equation ArPi{r,L) = 0. The exact form of the MFPT derived in Eq. (30 1 therefore provides 



straightforwardly an exact expression of the splitting probabilities. In particular, it admits a simple closed form in 
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150 200 



FIG. 7: Splitting probabili ty f or Levy flights on a network of size 200, for = 1.25, in discrete time (black line), compared to 
the theoretical expression ( 53 1 (red dashed line), and to the exact expression ( 52 1 where the first subleading term of the MFPT 



of Eq. (471 is taken into account (green dotted line). 
(-Pi,th) and the simulation (Pi,sim). 



The inset shows the difference between the two theoretical expressions 



the continuous limit, which reads 



Pi 



(-D 



l C(l-rf.„^) + C(l-c^■.,^)-C(l-rf^,f)-C(l-^^.,l-f) ^ 1 
2 2C(i-d^,^)-C(i-rf«„o)-c(i-c?^,i) ^ 



(53) 



For — 2 (Brownian case), we retrieve the classical expression, using (^{~n,x) = —Bn+i{x)/{n + 1), where _B„ are 
the Bernouilli polynomials {B2(x) — x"^ — x + 1/6), 



P,ix) 



i?2(f)+i?2(l-f)-i?2(i±^)-S2(i 



2 / 



2^2(0) +2^2(1)- 4^2 (i) 



= 1-x. 



In the generic case dw ^ 2, Eq. ( |53[ ) provides an exact expression for the splitting probability of a ID Levy process. 
We compare this exact solution in the continuous limit (53) with numerical simulations of Levy flights. The results 
are shown in Figs. [5] [6j and [7j indicating a rather slow convergence to the continuous limit as the system size 
grows. Clearly, the residual MFPT analyzed in previous sections yields a non vanishing correction to the splitting 
probabilities in the large system size limit. Fig. [7]shows that taking into account the residual MFPT in the evaluation 
of the MFPTs entering Eq. ( 52 ) , significantly improves the results for finite size L in the regimes of small r and small 
L — r. These results can be straightforwardly generalized to a larger number of targets using the formalism developed 
in Ref. 



MFPT FOR FRACTIONAL BROWNIAN MOTION 



Expressions (16) and (53) are valid for any continuous, compact, self-similar and Markovian ID random walk. 



Fractional Brownian motion (FBM) shares all these properties except for the Markov condition. FBM is defined as 
a continuous-time gaussian process with zero mean and stationary increments. More specifically, FBM is defined 
through the Langevin equation 



dx{t) 
dt 



m 



(54) 



for the position x{t), which is driven by the stationary, fractional Gaussian noise f (t) with (^(i)) = and long-ranged 
noise correlation [25| 
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FIG. 8: MFPT for an FBM, for several values of the Hurst exponent H = (0.40, 0.45, 0.50, 0.55, 0.60 and 0.65, from 

black [bottom] to light grey [top] as H grows). The continuous lines represents the simulations results, and the dashed lines 
the theoretical expressions following Eq. (16 1, where A is a fitting parameter. 



where we chose a unity diffusion constant. This process has the following covariance, 



(56) 



For subdiffusive processes (iJ < 1/2), the underlying fractional Gaussian noise (55) is anticorrelated, while it is 
positively correlated in the supcrdifFusive case 1/2 < < 1). 

To test numerically the applicability of our above results to non-Markovian processes, we performed simulations 
of FBM and compared the MFPT to a target for this process with expression (16). The results shown in Fig. [S] 
demonstrate that expression ( 16 1 with free parameter A provides a surprisingly good approximation of the MFPT 
for FBM, even though this process is highly non-Markovian (see also the discussion in Ref. [26 ). These numerical 
findings suggest that memory effects do not play an overly crucial role in the determination of MFPTs for FBM, and 
that the range of applications of the approach developed herein in practise extends to examples of non-Markovian 
processes. 



VI. CONCLUSION AND SUMMARY 

To conclude, we obtained a functional equation for the MFPT of a generic self-similar Markovian, continuous process 
to a target in a \D domain, and derived its exact solution. We showed that such a continuous limit of the MFPT 
is actually different from the large system size limit for discrete jump processes allowing leapovers. In the leapover 
case, the large system size limit of the MFPT involves non-vanishing corrections, which we call residual MFPT. This 
residual time can have important consequences in the context of both search processes and numerical simulations of 
first-passage times of random walks. We have investigated in detail the case of Levy flights, and validated our results 
by numerical simulations. We also demonstrated numerically, that our results apply with surprising accuracy to FBM, 
despite its non-Markovian nature. 
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